Operator splitting implicit integration factor methods for stiff reaction–diffusion–advection systems
نویسندگان
چکیده
For reaction–diffusion–advection equations, the stiffness from the reaction and diffusion terms often requires very restricted time step size, while the nonlinear advection term may lead to a sharp gradient in localized spatial regions. It is challenging to design numerical methods that can efficiently handle both difficulties. For reaction–diffusion systems with both stiff reaction and diffusion terms, implicit integration factor (IIF) method and its higher dimensional analog compact IIF (cIIF) serve as an efficient class of time-stepping methods, and their second order version is linearly unconditionally stable. For nonlinear hyperbolic equations, weighted essentially non-oscillatory (WENO) methods are a class of schemes with a uniformly high order of accuracy in smooth regions of the solution, which can also resolve the sharp gradient in an accurate and essentially non-oscillatory fashion. In this paper, we couple IIF/cIIF with WENO methods using the operator splitting approach to solve reaction–diffusion–advection equations. In particular, we apply the IIF/cIIF method to the stiff reaction and diffusion terms and the WENO method to the advection term in two different splitting sequences. Calculation of local truncation error and direct numerical simulations for both splitting approaches show the second order accuracy of the splitting method, and linear stability analysis and direct comparison with other approaches reveals excellent efficiency and stability properties. Applications of the splitting approach to two biological systems demonstrate that the overall method is accurate and efficient, and the splitting sequence consisting of two reaction–diffusion steps is more desirable than the one consisting of two advection steps, because CWC exhibits better accuracy and stability. 2011 Elsevier Inc. All rights reserved.
منابع مشابه
Operator splitting implicit integration factor methods for stiff reaction-diffusion-advection systems
For reaction-diffusion-advection equations, the stiffness from the reaction and diffusion terms often requires very restricted time step size, while the nonlinear advection term may lead to a sharp gradient in localized spatial regions. It is challenging to design numerical methods that can efficiently handle both difficulties. For reaction-diffusion systems with both stiff reaction and diffusi...
متن کاملKrylov single-step implicit integration factor WENO methods for advection-diffusion-reaction equations
Implicit integration factor (IIF) methods were developed in the literature for solving time-dependent stiff partial differential equations (PDEs). Recently, IIF methods are combined with weighted essentially non-oscillatory (WENO) schemes in [Jiang and Zhang, Journal of Computational Physics, 253 (2013) 368-388] to efficiently solve stiff nonlinear advection-diffusion-reaction equations. The me...
متن کاملEfficient and accurate splitting methods for time integration of multi-physics systems
Coupled multi-physics systems dominate computational models of advanced science and engineering applications. These systems give rise to challenges not typically faced by single-physics models: including interacting stiff and nonstiff processes, disparate time scales, increased nonlinearity and large data requirements. Historically, opposite approaches have been pursued by theoreticians and pra...
متن کاملNumerical Time Integration for Air Pollution Models
Due to the large number of chemical species and the three space dimensions, o -the-shelf sti ODE integrators are not feasible for the numerical time integration of sti systems of advection-diffusion-reaction equations @c @t +r (uc) = r (Kr c) +R (c) ; c = c(x; t); c 2 IRm; x 2 IR3 from the eld of air pollution modelling. This has led to the use of special time integration techniques. This paper...
متن کاملOn second derivative 3-stage Hermite--Birkhoff--Obrechkoff methods for stiff ODEs: A-stable up to order 10 with variable stepsize
Variable-step (VS) second derivative $k$-step $3$-stage Hermite--Birkhoff--Obrechkoff (HBO) methods of order $p=(k+3)$, denoted by HBO$(p)$ are constructed as a combination of linear $k$-step methods of order $(p-2)$ and a second derivative two-step diagonally implicit $3$-stage Hermite--Birkhoff method of order 5 (DIHB5) for solving stiff ordinary differential equations. The main reason for co...
متن کامل